A correlation between the equilibrium and transport properties 

of intercalation systems. 



E.V. Vakarin and J. P. Badiali 

Laboratoire de Electrochimie et Chimie Analytique 
ENSCP-UPMC, 11 rue P. et M. Curie, 75231 Cedex 05, Pans, France 

Employing the lattice gas model, combined with the linear elasticity theory, a cor- 
relation between the equilibrium and transport properties of intercalated species is in- 
vestigated. It is shown that the major features of the intercalation isotherms and the 
concentration dependence of the chemical diffusion coefficient can be well understood in 
terms of the change of the host volume in the course of intercalation. Theoretical predic- 
tions are compared to the experimental observations on PdH x , Li x WOs and Li-graphite 
systems. 
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I. INTRODUCTION 

Intercalation processes find their application in 
many technologically important domains, such 
as the design of hydrogen-storage systems [1], 
rechargeable high-energy batteries, electrochromic 
devices, (see Ref. [2] for a review), and supercon- 
ductors [3]. 

The insertion process can be viewed (at least, 
in some aspects) as an adsorption of guest parti- 
cles on the host lattice. For charged particles the 
ionic charge inside the matrix is compensated by 
the electrons. For that reason the electrochemical 
intercalation is also similar to a three-dimensional 
adsorption of neutral species. Based on this anal- 
ogy the intercalation is traditionally described 
within the lattice gas (LG) model. In this ap- 
proach all the properties (intercalation isotherm or 
capacity-concentration dependence) are connected 
with an ordering the guest on different adsorption 
sites of a rigid host lattice - the configurational 
transitions. 

However, the host response to the accommoda- 
tion of the guest species is not negligible. For in- 
stance, the insertion induces the stress into the 
host matrix. This may lead to segregation effects 
[4] or even to instabilities [5] of the host-guest sys- 
tem. Also a loading path is shown [6] to influ- 
ence the guest uptake efficiency. Due to this the 
host may undergo an expansion or local distor- 
tion. Typical examples are the hydrogen sorption 
by metals [6-9] , or the intercalation of Li ions into 
layered materials [10-13]. Quite often the host un- 
dergoes structural transformations [14-17] due to 
intercalation. These are the structural transitions. 
In such cases the standard LG approach also oper- 
ates with several sublattices [10,17] corresponding 
to each configuration of the host. However, such 
a restructuring suggests that the elastic effects 
should be taken into account. A general thermo- 
mechanical theory of the stress-composition in- 
teraction is developed by Larche and Cahn [18]. 
This implies the existence of a coupling between 



the elastic properties of the host material and the 
structure and/or the dynamics associated with the 
guest species. In other words the configurational 
and structural transitions should be considered to 
be coupled. In our previous papers we have inves- 
tigated such a coupling for two-dimensional [19] as 
well as for three-dimensional systems [20-22] . The 
main advantage of our approach is its capability 
in describing various (microscopically different) in- 
tercalation systems within a common framework. 
Despite of the fact that some microscopic details 
are omitted, the theory agrees well with exper- 
imental data on layered Li x TiS2 and crystalline 
Li x WOs and Na x WOs compounds at equilibrium 
conditions. 

In this study we focus on a correlation between 
the equilibrium properties (isotherms and differ- 
ential capacity) and the concentration dependence 
of the chemical diffusion coefficient. Our purpose 
is to determine the impact of the thermodynam- 
ics to the kinetic properties. Following our previ- 
ous works we combine the LG model, describing 
the configurational transitions and the linear elas- 
ticity theory that accounts for the loading mech- 
anism and the volume dilatation in the course of 
the structural transitions (expansion or restructur- 
ing). This allows us to derive an effective chemical 
potential, involving the concentration dependent 
stress and strain fields. Then, assuming that the 
diffusion flux is proportional to the gradient of the 
chemical potential, we investigate the behavior of 
the chemical diffusion coefficient. 



II. MODEL 

A. Microscopic formulation 

The host material is described as a three- 
dimensional lattice of adsorbing sites with their 
positions given by a set vectors (rj. Due to 
the elastic properties of the real host each site 
(for instance, an interstitial site) of this "auxil- 



iary" lattice may deviate from its equilibrium po- 
sition r°, such that we deal with the displacements 
Ui = Ti — r?. Therefore, the host properties are de- 
scribed by the Hamiltonian i?#[{uj}]. Note how- 
ever that a connection between the elastic proper- 
ties of the real matrix and those of the adsorbing 
lattice, is not straightforward [20,21]. 

A distribution of the intercalants on the host 
sites is given by a set of occupation numbers {ti}, 
with ti = or ti = 1. The guest subsystem is char- 
acterized by the chemical potential \i and the near- 
est neighbor interaction parameter W (at equilib- 
rium positions). For the electrochemical insertion, 
chemical potential /i gives a deviation of the elec- 
trode potential — eV from its standard value E . 
In this way an arrangement of the intercalants on 
a rigid host lattice is governed by the LG Hamil- 
tonian 

H G = W^2t i t j -vJ2 ti w 

ij i 

that describes the configurational transitions of the 
intercalated species. These could be the droplet 
formation for attractive interactions (W < 0) or 
the order-disorder transition for repulsive interac- 
tions (W > 0). In the latter case one introduces 
the sublattice concentrations as appropriate for the 
symmetry of a given system. 

The coupling between the host and the guest is 
given by Hamiltonian Hc[{ui}, {ti}} which takes 
into account a dependence of the binding energy 
on the site displacement and also the pairwise in- 
teraction between the guest particles through the 
host lattice. The overall Hamiltonian is now writ- 
ten as 

H = H H [{ui}} + H G [{ti}} + H c [{ui}, {*<}] (2) 

The free energy F corresponding to the above 
Hamiltonian is given by 

F = F H +F G {x)+F c {x) (3) 

where x is the intercalant concentration. Here Fh 
is the host free energy in the absence of intercala- 
tion, F G is the guest free energy in a case of the 
rigid host lattice. The latter term can be calcu- 
lated using the mean field approximation to obtain 
the well-known relation [2,10] for the chemical po- 
tential 

f io(x) = -eV-E = cWx+^ln(^—^ (4) 

where c is the coordination number of the host 
lattice and (3 = 1/kT. The coupling term j3F G = 
— ln(((e _/3ffc )( u .))( t .)) requires the averaging over 
the displacements and the occupation numbers, 
calculated with the reference terms Hu and Hq. 
In fact, this is an infinite series including the 



correlations of all orders in the reference state. 
The main problem is to specify Hc[{u.i}, {U}} 
coherently with the host symmetry and elastic 
properties. It is known that real host materi- 
als have rather complicated elastic properties (e.g. 
a strong anisotropy). Therefore, only some sim- 
plified model calculations are expected to give 
tractable results. On the other hand, such predic- 
tions (e.g. a rigid plane model [23]) do not exhibit 
quantitative agreement with experimental data. 

B. Approximation 

The perturbative scheme developed previously 
[20,21] has shown that the coupling term is con- 
cerned with a concentration dependence of the host 
response to the intercalation. This involves at least 
two effects. First is a renormalization [2,24] of the 
net pair interaction between the intercalants due to 
their indirect interaction through the matrix. The 
second is a change of the host volume upon inser- 
tion of the guest species. Depending on the host 
nature, a stress field may result if the lattice is not 
totally free to relax. In this situation it seems rea- 
sonable to estimate the coupling term based on the 
continuum theory of elasticity with the concentra- 
tion dependent stress and strain fields. Then the 
host-guest free energy is a sum of the lattice gas 
and elastic part 

F(x) = F LG (x) + F el (x) (5) 

where Flg{ x ) is the configurational (lattice gas) 
part, in which the pairwise interaction is renor- 
malized [2] due to the interaction through the ma- 
trix. This gives the chemical potential /xo(x) (see 
eq. (4)) with a new interaction constant W. The 
elastic part is approximated by the free energy of 
a strained isotropic body under a loading stress 
o(x). Since the strain is assumed to be purely di- 
latational, we operate with traces e and a of the 
corresponding tensors. 

F el (x) = ^e(x) 2 -a{x)e{x) (6) 

with A being the effective elastic constant, inde- 
pendent of the concentration. The total stress 
S = S(x) is given by 

dF i 

S(x) = ^=As(x)-a(x) (7) 

Therefore, we have two stress contributions. The 
internal, or self-stress Ae(x) corresponds to the 
host reaction to the guest insertion. The second 
term a(x) describes a loading procedure, that may 
include the sample clamping or other effects which 
are not directly related to the strain. It is impor- 
tant that, in general, a(x) is a function of x (not 



a function of s(x)). For instance, if the sample 
is clamped such that e{x) = 0, then we have a 
stress accumulation proportional to the concentra- 
tion a(x) oc x. 



III. EQUILIBRIUM PROPERTIES 

The guest chemical potential ji{x) is given by the 
concentration derivative of the total free energy. 



H{x) = (i (x) + S(x) 



de(x) dcr(x) 



dx 



dx 



e(x) (8) 



Here the second term involves the so-called chem- 
ical expansion coefficient de(x)/dx 1 while the last 
term is associated with the loading path. It is seen 
that the intercalation level depends on an inter- 
play of the internal stress and the loading stress. 
The latter could be small, but its concentration 
derivative is not necessary small, so that the load- 
ing path may induce serious consequences [6]. In 
particular, for a given material (A) and a suitable 
loading path a(x), there may be a cancellation of 
the last two terms in eq. (8) in a given domain 
of x. This explains why in some cases the purely 
configurational description {u{x) = ^o( x )) works 
well. 

For comparison to other theoretical approaches 
it is convenient to express the strain function in 
terms of other relevant quantities. The strain can 
be measured as a volume dilatation or as a change 
of the interlayer spacing in the course of insertion. 
If V{x) is a composition dependent sample volume 
(or interlayer spacing for layered compounds), then 
by definition 



e{x) = 



V(x) - V(0) 



Sp(x) 



(9) 



where <5 = [V(l) - V(0)]/V(0) is the relative vol- 
ume variation and p(x) = [V(x) — V(0)]/[V(1) — 
V(0)] is a modulating function. The latter varies 
in between and 1. If p{x) is a linear function, 
then the system is said to obey the Vegard's law 
[25]. 

If the host-guest system forms a solid solution, 
then the guest partial molar volume V m (x) is re- 
lated to the total volume 



V(x)=xV m (x) + V(0), 



(10) 



where V(0) is the initial host volume (at x = 0). 
Then we can find the following relation for the 
strain, the volume variation and partial molar vol- 
ume 



V(0) 



V(0) 



(11) 



If the loading is composition independent <r{x) = a 
and the sample volume changes linearly (V m = 



const, e{x) = xV m /V(0), like in the PdH x a- 
phase), then we recover the well-known result [9] 



H(x) = u (x) 



"-"m)m) (12) 



The term liner in the concentration can be viewed 
as an additional interaction which can be com- 
bined with cWx in /io(x). Although this additional 
term is proportional to the elastic constant A, it is 
physically different from the interaction of elastic 
dipoles [2,24] which is already absorbed into W 
(see the discussion after equation (5)). This inter- 
action takes place even if the sample is perfectly 
clamped (e(x) — 0), while our term represents a 
cooperative effect due to e(x) ^ 0. 

In order to analyze the role of elastic effects in 
the phase behavior of the guest species it is in- 
structive to consider a simple example of vanish- 
ing loading stress a(x) = 0, S(x) = Ae(x). In 
this case there is no difference between the volume 
expansion or contraction. Then we arrive at 



u(x) = n (x) + Ap(x) 



dp(x) 
dx 



(13) 



where A = AS 2 . In addition, p(x) is taken to be 
that in the so-called layer rigidity model [25] which 
describes a deviation from the Vegard's law in lay- 
ered compounds 



p{x) = 1 - (1 - xf 



(14) 



where q is related to the rigidity of the host lay- 
ers. The first order phase transition (of liquid-gas 
type) is manifested by a singularity in the differ- 
ential capacity C — dx/d([3u) 



fiC 



d Vo(x) , t ( dp{x) 



dx 



A 



\ dx 



A 



d 2 p(x) 
dx 2 



(15) 



It is clear that for q < 1 both derivatives of p(x) 
are positive and the only possibility for the criti- 
cal behavior is the case of attractive interactions 
W < 0, i.e. when the configurational part 



du (x) _ 1 + f3cWx(l - x) 



dx 



x(l — x) 



(16) 



becomes negative. Then the last two terms in 
the denominator of (15) rescale the critical tem- 
perature and shift the critical composition from 
x = 1/2. For q > 1 the term d 2 p/dx 2 becomes 
negative and the criticality may appear due to the 
elastic effects even if W = 0. The coexistence 
curves for this case, obtained by a numerical anal- 
ysis of the divergence in C, are plotted in Fig. 1. 



It is seen that the critical temperature (the maxi- 
mum) increases with increasing q. Simultaneously 
the critical concentration decreases. The elastic ef- 
fects break the hole-particle symmetry of the prob- 
lem and consequently the curves are not symmetric 
and are not centered around x = 1/2. 

In our previous works [20-22] we have ana- 
lyzed another extreme case of vanishing total stress 
S (x) — > and the loading stress proportional to 
the concentration a(x) = gqx. Then 



(17) 



where 7 = Sa - Note, however that, for A = const, 
the above result is valid only for a linear behavior 
of p(x) = px, while is an approximation for an ar- 
bitrary p(x) . In reality the elastic constants could 
depend on the concentration, for instance the bulk 
modulus of Pd is reduced [8] up to 20% due to the 
hydrogen sorption. Thus one can easily imagine 
a situation when A(x)s(x) « &(x) and then the 
approximation (17) is applicable for a nonlinear 
p(x). In general, for a non-Vegard's behavior and 
the linear loading a(x) = a x we deal with 



ji(x) = no(x) + [Ap(x) 



dp{a 



dx 



!P{x) (18) 



The strain can be measured as a volume dilata- 
tion or as a change of the interlayer spacing during 
the transitions between different phases (staging 
in graphite, restructuring in Li x WOz, a — (3 tran- 
sition in PdH Xl etc). In any case the lattice pa- 
rameters do not obey the linear Vegard's law [25], 
exhibiting a sharp change near the transition com- 
positions x® n . This can be described by the follow- 
ing approximation [20-22]. 



P{x) = \ 



l + ^Pn tanh[a„(a; 



(19) 



where n counts the number of phase boundaries. 
In LixWOz, n = 1,2 corresponds to monoclinic- 
tetragonal and tetragonal-cubic transitions, re- 
spectively. For Lz-graphite systems, n corresponds 
to the stage transfer boundaries [12,13]. The set of 
rigidity parameters a n controls a local slope and 
p n are the weights corresponding to each phase, 
such that J2 n P n = 1- This * s consistent with 
experimental observations [14] indicating that the 
structures are not completely pure, but contain 
some features that indicate a mixing of phases. 
The present form of p(x) corresponds to contin- 
uous structural transitions, but can be easily mod- 
ified to take into account the jump-like behavior 
(like for the staging [2]). Nevertheless, in that case 
the concentration derivatives of p(x) would be sin- 
gular, inducing the singularities in the thermody- 
namic quantities (as it should be in the neighbor- 
hood of a phase transition). Note however that the 



experimental dependencies are usually smoothed 
due to a finite concentration resolution. Then it is 
difficult to distinguish between the step-wise vari- 
ation and a sharp (but continuous) transition. For 
that reason the criticality criteria, determined on 
the ground of rigorous statistical mechanical argu- 
ments (u(x) loops or a divergent capacity C(x)), 
only approximately conform to experimental data. 
Therefore we stay with the continuous p{x), per- 
forming the fitting to the experiment. 

Following this methodology, in Fig. 3 we display 
the intercalation isotherm (inset) and the capacity 
curve, comparing them to the experimental data 
[14] for Li x W0 3 . The fitting (eq. (18)) is per- 
formed assuming that the effective pair interaction 
inside the matrix is repulsive for any x. The curve 
modulation results from the strain behavior p(x), 
associated with the volume dilatation. As is dis- 
cussed previously [21,22], the peculiarities (inflec- 
tion points and the peaks in C(x) = dx/d{{3u)) 
mark the boundary between different host sym- 
metries. From the experimental point of view, 
the peak height is usually associated with a tran- 
sition sharpness. Nevertheless, as we have dis- 
cussed above, the critical behavior of the model 
itself should be studied separately. 



IV. CHEMICAL DIFFUSION COEFFICIENT 

In many cases there are two mobile species in- 
side host matrices: the host electrons and the guest 
particles. Therefore, the transport is characterized 
by their flux densities J e and 3 G , related to the 
gradients of the corresponding chemical (or electro- 
chemical) potentials u e and a G . 



J e = L ee Vu e + L eG V/i G 



J G = L GG V/i G + L Ge \7u e 



(20) 



(21) 



where L a \, arc phenomcnological transport coeffi- 
cients, which are scalar quantities if the host is 
spatially isotropic. In matrices with a metallic 
conductivity the electrons are much more mobile 
than the guest species. Due to this the guest ionic 
charge is compensated by the electrons. For the 
same reason the electrons reach a uniform (equi- 
librium) distribution much faster than the inter- 
calants ( V/i e = on the time scale when V/i G ^ 
0). In addition, usually L GG >> L eG [2]. There- 
fore one may focus on the guest flux 



J = -M{x)xVu = -D(x)Vx 



(22) 



where M(x) is the mobility coefficient that must 
contain a blocking factor M(x) = D (l — x), deter- 
mining the concentration-dependent chemical (or 
collective) diffusion coefficient D(x). A compre- 
hensive review of theoretical approaches to the 



determination of the chemical diffusion coefficient 
can be found in Ref. [26]. In the simplest approx- 
imation we deal with 



D(x) = D x(l - x) 



d/3[i 
dx 



D x(l-x)/C(x) (23) 



Under the assumptions above, the intercalant 
transport is described by the diffusion equation 
with an effective diffusion coefficient D(x). The 
latter requires an information on the intercalation 
isotherm n(x). Then, based on the equilibrium 
properties, at a given concentration gradient, one, 
at least in principle, can solve the kinetic problem. 

Note, that such a simple scheme should be mod- 
ified when the host conductivity changes remark- 
ably upon the intercalation. Then the electronic 
impact should be taken into account. Also, we 
do not consider other driving forces, like external 
field or stress gradients, assuming that the guest 
concentration, x, is the only independent variable. 
And finally, near the phase coexistence (e.g, stag- 
ing) a system becomes non-uniform because of the 
phase boundaries. Then the formulation of the ki- 
netic problem must be coherent with the theory of 
critical phenomena (see [30] for a recent review). 

Starting from (8) and (23) we obtain 



D{x) = D (x) + D el {x) 



(24) 



where the first term is the standard LG part 
Do(x) = D x(l — x)^r and the elastic contribu- 
tion is given by 



D P 



B(x) 



S— + a(—\ 

dx 2 \dx J 



^ de da 
dx dx 



fa 

dx 2 ' 



with B{x) = Dqx(1 — x). The diffusion coeffi- 
cient involves several competing factors (strains, 
stresses and their concentration derivatives). In 
general it is not trivial to see whether D increases 
or decreases with x. The situation is even more 
complicated for systems in which the lattice spac- 
ing docs not obey the linear Vegard's law, but has 
inflection points separating different phases (e.g. 
staging in graphite or restructuring in Li x WO^). 
Then the derivatives above may change sign with 
the concentration. Also we see again that the load- 
ing path has a significant contribution. 

For PdH x (a-phase) we recover the well known 
result [8,9]. The concentration induced internal 
stress Sq(x) = AxV m /Vo increases the diffusion 



D = D (x)+B(x)A[^ 



(26) 



where D (x) is the stress free contribution (it cor- 
responds to /xo(x), i.e. the lattice gas description, 
see above). Note that the non-local stress effects 
[7,8] are not discussed here. 



Since, according to (23), the diffusion coefficient 
is just an inverse of C(x), D as a function of x 
has minima (Fig. 3), corresponding to the peaks 
in the capacity curve. For WO3 (Fig. 3) the the- 
ory works well near the first minimum, but overes- 
timates the diffusion at higher x. Assuming that 
eq.(23) is valid, we have also inverted the capacity 
data (Fig. 2) in order to demonstrate that there is 
no "hidden" errors in the fit for C(x). Since the 
equilibrium characteristics (Fig. 2) are predicted 
with a reasonable accuracy, then one of the source 
for the discrepancy could be our simplified assump- 
tions on the mobility. For instance, one might sup- 
pose that eq.(23) does not work for higher x be- 
cause the hoping D depends on the elastic prop- 
erties. By the analogy with the two-dimensional 
diffusion [26], Do is related to the lattice spacing. 
The latter decreases with x for WO3 [15], so that 
the resulting D(x) would also decrease. Never- 
theless the volume change does not exceed several 
percents [2,15]. This cannot explain the two or- 
ders of magnitude difference in Fig. 3. The lattice 
anisotropy seems also to be irrelevant because the 
final (high x) structure is cubic, so the isotropic 
approximation for the diffusion and elastic proper- 
ties is reasonable. 

It is known [27] that Li x WOz changes its elec- 
tronic conductivity near the structural instability 
composition. Then the ionic transport must cor- 
relate with the transport of the neutralizing elec- 
trons. Similar effects occur in superconductors [3]. 
Therefore, the electronic mobility impact (that is 
absent from our approach) to the current should 
/„be taken into account by including the additional 
^ driving force V^ e = dfi e /dx\7x. This contribution 
can be included into our approach, assuming some 
model for the electronic structure, e.g. - the rigid 
band model. This would correct our predictions 
for small x (according to [27], for x < 0.1), but our 
high- a; estimation would remain unchanged. We 
believe that the electronic effects are not responsi- 
ble for the decrease of D(x) at high concentrations. 
Based on the fact [27] that the electrical properties 
of Li x WOs become more and more metallic with 
increasing Li content, we may expect that our ap- 
proximation V/x e = becomes more reliable with 
increasing x. Moreover, if the conductivity effects 
were important, then the ionic diffusion would in- 
crease like this occurs in amorphous WO3 films 
[28]. 

On our opinion we deal with an interplay of sev- 
eral facts. Although the model works well in pre- 
dicting the equilibrium properties, nevertheless it 
is probably too simple for the description of the 
kinetics. In fact, the transport is assumed to take 
place in an infinite lattice through the hoping of 
neutral species (Li + + e), ignoring the kinetics 
other relevant phenomena, like the charge transfer 
(metal/host, guest/host), the exclusion (permse- 



lectivity) effect [20], the formation of the passive 
layer at the host/electrolyte boundary, etc. On the 
other hand, the experimental data on the diffusion 
are extracted from direct measurements employing 
a theoretical model, which could differ (in some 
details) from the one introduced here. Concern- 
ing the order of magnitude estimations, it should 
be noted that the experimental data for the same 
substance are usually differ significantly depend- 
ing on the sample preparation, its size and the 
scanning time interval. Therefore, the quantita- 
tive description of the ionic transport is a delicate 
problem requiring mutual theoretical and experi- 
mental efforts. Nevertheless, our estimation of the 
elastic effects is qualitatively correct in predicting 
that the system characteristics changes coherently 
with the strain p(x), and therefore, positions of the 
minima (Fig. 3) x^ are predicted correctly, imply- 
ing a correlation between the D{x) behavior and 
the strain development. 

Similar situation takes place for iz-graphitc. 
In Fig. 4 the normalized interlayer spacing for 
graphite in the presence of intercalated Li is plot- 
ted MS iX function of the guest composition. The 
fitting is preformed using eq. (19) under a suitable 
choice of the parameters p n , a n , x^- The experi- 
mental data are taken from Refs. [12,13]. Our fit 
is rather reasonable, "catching" the fact that the 
stage transfer compositions x n approximately cor- 
respond to the inflection points in p(x), where the 
average interlayer spacing changes sharply. Note 
that we do not discuss the fine structure of the 
phase diagram [12], such as dilute and liquid-like 
phases for the same stage. These features are re- 
lated to the the guest in-plane ordering which is 
indistinguishable in the behavior of the interlayer 
spacing. 

Having a reliable approximation for p{x) we cal- 
culate the diffusion coefficient in the framework 
outlined above. The fitting to the experimental 
data [31] is shown in Fig. 5. It is seen that the 
diffusion slows down significantly near the stage 
transfer compositions Although the magni- 

tude of D{x) is close to the experimental results, 
the width of the minima is underestimated. Note 
that the fitting is not optimized, that is, we did not 
try to find an optimal set of the parameters, which 
gives equally good agreement for the isotherm ^i(x) 
and D(x). Nevertheless, as in the case of Li x WOs 
(Fig. 3), there is a clear correlation between the 
diffusivity and the strain behavior. 



V. CONCLUSION 

In summary, our approach implies that the main 
equilibrium and transport features of the intercala- 
tion systems (which differ in their microscopic de- 
tails) can be well understood in terms of a concen- 



tration dependence of the hydrostatic parts of the 
stress and strain fields, associated with the internal 
and loading effects [6] . It is shown that the elastic 
effects may induce the critical behavior even if the 
"direct" interaction between the guest species is 
absent. Although the ionic transport is supposed 
to be equivalent to the diffusion of neutral (ion plus 
electron) species, the approach is flexible enough 
to incorporate other transport mechanism (migra- 
tion, or electronic mobility effects). 

The theory gives a quantitative description of 
different insertion processes, involving the volume 
expansion {Li x TiSi (see [20]) , a — PdH x ), staging 
(Li-graphite) or restructuring (Li x WO^). For all 
these processes the theory implies a correlation be- 
tween the intercalation isotherm and a concentra- 
tion dependence of the diffusion coefficient. The 
latter exhibit a set of characteristic minima, re- 
lated to the boundaries between different phases 
(like different symmetry phases of Li x W03, a — (3 
transition in PdH Xl etc). However the experimen- 
tal dependencies are usually smoothed due to a 
finite concentration resolution. For that reason 
the criticality criteria, determined on the ground 
of rigorous statistical mechanical arguments (e.g. 
D{x) — at the transition concentrations) only 
approximately conform to the experimental data, 
exhibiting a sharp (but finite) decrease of D. 

For a non-Vegard's strain variation the diffusion 
coefficient is a nonlinear (and nonmonotonic) func- 
tion of the concentration (25). Therefore the dif- 
fusion equation for the concentration profile would 
be strongly nonlinear. Then one can expect a 
rather complicated space-time variation, including 
, for instance, oscillations [6] and other nonlinear 
effects. 

Our results may have implication in various do- 
mains related to the insertion process, like hydro- 
gen sorption [1], electrochemical intercalation, im- 
purities in alloys, layered superconductors [3], vol- 
ume transitions in hydrated gels [32], etc. 
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FIG. 1. Liquid-gas coexistence curves for the layer 
rigidity model in the absence of the direct interaction 
W = 0. 



FIG. 2. Differential capacity and the voltage (inset) 
for crystalline Li x WOs- The symbols correspond to 
the experimental data [12]. The parameters pi — 0.7, 
p 2 = 0.3, qi = q 2 = 15, x\ = 0.05, x° 2 = 0.27, 
/?7 = 1.8, /3A = 0.01, j3W = 2.2. 



FIG. 3. The chemical diffusion coefficient for crys- 
talline Li x WOz. The symbols (up triangles) corre- 
spond to the experimental data [24]. The parameters 
are the same as for the previous figure. 



FIG. 4. The average interlayer spacing for 
Li-graphite. The experimental data are extracted from 
[12,13] . The theoretical fitting is performed using equa- 
tion (19), where pi = 0.28, ati = 30, x\ = 0.04, 
p 2 = 0.22, a 2 = 20, x% = 0.25, p s = 0.5, q 3 = 10, 
xl = 0.75 



FIG. 5. The chemical diffusion coefficient for 
Li-graphite. The experimental data are extracted from 
[31]. The fitting is done using equations (18) and (23) 
where (3cW = -0.1, 7 = 0.5 A = 0.1 
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